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Decay of an isolated monopole into a Dirac monopole configuration 
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We study numerically the detailed structure and decay dynamics of isolated monopoles in condi¬ 
tions similar to those of their recent experimental discovery. We find that the core of a monopole 
in the polar phase of a spin-1 Bose-Einstein condensate contains a small half-quantum vortex ring. 
Well after the creation of the monopole, we observe a dynamical quantum phase transition that 
destroys the polar phase. Strikingly, the resulting ferromagnetic order parameter exhibits a Dirac 
monopole in its synthetic magnetic field. 


Introduction —The significant roles played by topolog¬ 
ical defects in nature and their appearance in various 
physical contexts mm have sparked numerous theoreti¬ 
cal and experimental studies. The precise control over ex¬ 
perimental parameters and the ability to image the quan¬ 
tum mechanical order parameter directly render Bose- 
Einstein condensates (BECs) a unique platform to create 
and observe different types of topological defects. In par¬ 
ticular, BECs with spin degrees of freedom may host a 
rich variety of defects due to many possible order param¬ 
eter manifolds and symmetries [sHini. In these systems, 
topological defects can either be created in a determinis¬ 
tic manner using precisely controlled magnetic and laser 
fields pTlUb] , or they can form spontaneously, for exam¬ 
ple when the condensate is rapidly quenched a through 
quantum phase transition [min]. The experimentally 
realized topological structures in BECs to date include 
singly and multiply quantized vortices [n na mHin], 
solitons and vortex rings [H], skyrmions [m |22], po- 
lar core vortices m, coreless vortices [23l [24], vortex- 
antivortex superpositions [25], solitonic vortices [26] and 
monopoles [ISIEI]. 

The Dirac monopole configuration created in Ref. m 
is an analogue of the classical magnetic point charge 
considered by Dirac in the context of quantum mechan¬ 
ics [28]. It manifests itself as a point-like singularity in 
the so-called synthetic magnetic field [29] , which is an ef¬ 
fective gauge field for the scalar part of the order param¬ 
eter arising naturally from its spin degrees of freedom. 
In agreement with Dirac’s original work [28], this kind of 
monopole induces in the condensate order parameter a 
nodal vortex line with vanishing particle density extend¬ 
ing from the location of the monopole to the boundary 
of the atom cloud. Thus the ferromagnetic order param¬ 
eter supporting the Dirac monopole is energetically and 
dynamically reminiscent of a line defect. Critically, there 
is no topological point defect in the order parameter it¬ 
self, as it is in a configuration topologically equivalent to 
the ground state. Indeed, the second homotopy group [2] 
for the ferromagnetic order parameter space is trivial and 
topological point defects are not permitted. Point defects 
may exist in the polar phase of a spin-1 condensate, how¬ 


ever, as the second homotopy group for the polar order 
parameter space, Gp = [S‘^ x t/(l)]/Z 2 [3Ql l3T] , is iso¬ 
morphic to the additive group of integers, 7r2(Gp) = Z. 

Topological point defects in the polar phase of a 
^^Rb spin-1 BEC have been recently realized experimen¬ 
tally m- No nodal lines or other physical line-like ob¬ 
jects are attached to this monopole, and we therefore 
refer to it as an isolated monopole. As Ref. EH focused 
only on the first observation of the isolated monopole 
within a finite experimental resolution, there has been 
no detailed study of the fine structure of the defect. Eur- 
thermore, both the isolated monopole [5] and indeed the 
entire polar phase of the ^^Rb condensate [17] are ex¬ 
pected to be unstable at low magnetic fields, prompting 
a study on the evolution of the isolated monopole after 
its creation. 

In this Letter, we present computational results of the 
fine-grained structure and decay dynamics of the iso¬ 
lated monopole. We observe that the core of the cre¬ 
ated monopole contains a small half-quantum vortex ring 
as expected based on energetic arguments studied in 
Ref. [5]. We show that beyond the experimentally ac¬ 
cessed time scales EH the polar order parameter evolves 
into a ferromagnetic order parameter, accompanied by 
the decay of the isolated monopole into a Dirac monopole 
in the resulting synthetic magnetic field. Importantly, 
this decay arises naturally from the thoroughly under¬ 
stood physics of the atom cloud without any phenomeno¬ 
logical damping terms. Quantitatively matching dynam¬ 
ics are observed for both ferromagnetic and antiferromag¬ 
netic spin-spin interactions. 

Theoretical background —The mean-field order param¬ 
eter of a spin-1 condensate can be expressed as '0(r) = 
Y^n(r)^(r), where n(r) is the particle density and ^(r) 
is a three-component complex-valued spinor such that 
^(r)’l'^(r) = 1. The temporal evolution of the order pa¬ 
rameter at low temperatures is accurately described by 
the differential equation 

= {/i(r) + n(r)[co (1) 

+ C2S(r) ■ F] - irn^(r)}V’(r), 
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where h{r) is the single-particle Hamiltonian, F = 
{Fx^Fy^Fz) is a vector composed of the dimensionless 
spin-1 matrices, F is the three-body recombination rate, 
and S(r) = ^(r)’t’F^(r) is the local average spin. The 
coupling constants characterizing the atom-atom inter¬ 
actions are given by cq = 47r/i^(ao + 2a2)/(3m) and 
C 2 = 47rfi^(a2 — ao)/(3m), where a/ is the 5- wave scat¬ 
tering length corresponding to the scattering channel 
with total two-atom hyperfine spin /. The single-particle 
Hamiltonian is given by 

h(r) = - / (2m) FV{y) - ji (2) 

+ 5FMBB(r, t)-F + g[B(r, t) ■ F]^, 


where m is the mass of the atoms, H(r) is an external 
optical trapping potential, /i is the chemical potential, gp 
is the Lande ^-factor, is the Bohr magneton, B is an 
externally applied magnetic field, and q is the strength of 
the quadratic Zeeman shift [32]. We assume that H(r) = 
m[co’^(x^ +co’^z^]/2, where ujr and ooz are the radial 

and axial trapping frequencies, respectively. 

In the pure polar phase with S(r) = 0, the order pa¬ 
rameter can be expressed in the basis of the 2 ;-quantized 
spin states {|1) , |0), | —1)} as [33ll34] 


V’(r) = 


Vnjr) 


Jfpir) 


V2 


-ax[r)+tay[L) 

V^d:,{r) 

d^{r) + idy(v) 
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Thus in the Cartesian basis the polar order parameter 
reads '0(r) = where d = {dx^dy^dz)^ is 

a real unit vector known as the nematic vector. In the 
case of non-zero average spin, we investigate the nematic 
order through the magnetic quadrupole tensor [33] 


Qab — 




( 4 ) 


where are the components of the spinor in the Carte¬ 
sian basis. For (S(r)) ^ 0, the vector d is defined as the 
eigenvector corresponding to the largest eigenvalue of Q. 

Methods —The nematic vector behaves identically to 
the average spin under rotations in spin space, and hence 
it also follows adiabatic changes in the external magnetic 
field. Consequently, the method originally developed in 
Ref. |8] for the adiabatic creation of Dirac monopoles in 
the ferromagnetic phase can be used to create isolated 
monopoles in the polar phase, as realized in Ref. [27|. In 
brief, the condensate is subjected to an external magnetic 
field B(r, t) = Bq(r)+Bb(t), where Bq(r) = 6q(xx+^y- 
22 ;z) is a quadrupolar magnetic field with gradient and 
Bb(t) = Bz{t)z is a spatially homogeneous bias field. In 
the beginning of the simulation, the condensate is in the 
spin state |0), yielding a nematic vector d(r) = z. At the 
initial bias field Bz = 1 G, the field gradient is linearly 
ramped from zero to = 3.7 G/cm in 10 ms and the 
bias field is subsequently decreased to Bz = 10 mC in 
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FIG. I. (a) Horizontally (b) and vertically integrated parti¬ 
cle densities for the indicated waiting times after the creation 
ramp which produces the the isolated monopole. Different 
colors correspond to particles in different ^-quantized spin 
states with the color and intensity scales given in the bottom¬ 
most panel. The peak column density is hp = 2.7 x 10^^ cm“^ 
and the field of view is 15.5 x 15.5 giM in each panel. The 
white arrows indicate the location of a vortex line shown more 
clearly in Fig. 
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10 ms. The monopole is created by decreasing the value 
of the bias field linearly to zero at a rate B^it) = —0.25 
G/s. This part of the control protocol is referred to as the 
creation ramp and it ideally results in d = Bq(r). After 
the creation ramp the temporal evolution continues with 
the quadrupole field and optical trap intact. 
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FIG. 2. (a) Spin density and (b) particle density 

of the condensate 200 ms after the creation ramp. The shown 
density range is [nmin, nmax] = [1.0,4.8] X 10“^ N/al. The 
data in both panels is shown only for |z| < 1.5ar. The spin 
density is well depleted along the vortex, whereas the particle 
density is only partially depleted. 


The initial particle number is 2.1 x 10^ and the optical 
trapping frequencies are uJr = 27r x 124 Hz and uJz = 
27r X 164 Hz. We take the other parameters according 
to ^^Rb such as the literature values for the atom loss 
arising from the three-body recombination T = h x 2.9 x 
10“^^ cm^/s [35l |36], the quadratic Zeeman shift q = 
27rh X 2.78 MHz/T [T7], and the scattering lengths ao = 
5.387 nm and 02 = 5.313 nm. The computational volume 
considered is 24x24x24a^, where {mujr) = 

1 /im and the corresponding size of the computational 
grid is 200x200x200 points. In order to enhance the 
numerical emulation of the experimental conditions EH, 
we add spatially uncorrelated complex noise to the spinor 
components at each grid point prior to the creation ramp. 
The amplitude of the noise is uniformly distributed to 
introduce 0-1% fluctuations in the local particle density 
for each spinor component. 

Results —We numerically integrate Eq. § and apply 
the control protocol described above with the initial con¬ 
dition d = z. Figure shows the resulting spin-contrast 
images of the condensate particle densities for different 
waiting times after the creation ramp. As in Ref. 
the condensate particle densities just after the creation 
ramp are in good agreement with Eq. and d = Bq(r). 
Thus the particles almost entirely reside in the so-called 
neutral state which corresponds to the zero eigenvalue of 
the local Zeeman Hamiltonian. However, 50 ms after the 
creation ramp, the polar phase has noticeably decayed 
towards the local ferromagnetic strong-held seeking state 
(SFSS), i.e., the spin state that minimizes the local Zee- 
man energy. The ferromagnetic phase is hrst visible at 
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FIG. 3. (color online) Projection of (a-c, e, f) the nematic vec¬ 
tor d (t) and (d, g) average spin S (t) of the condensate onto 
(a-d) the xz- and (e-g) xy-planes for different indicated wait¬ 
ing times after the creation ramp. The background colormap 
shows (a-d) the y component, (e,f) the azimuthal component, 
and (g) the z component of the spin. The core of the Alice 
ring is shown with black dots in panel (a) which shows a mag- 
nihcation of the center region of panel (b). The white arrow 
indicates the location of a vortex line shown more clearly in 
Fig. [2] 


the top and bottom edges of the condensate and extends 
gradually until the condensate resides almost entirely in 
the SFSS. Qualitatively similar results are obtained in 
simulations without the added noise (data not shown). 

It is well known that the SFSS in the quadrupole field 
Bq(r) corresponds to the synthetic magnetic field of a 
Dirac monopole [8]. However, the final order parame¬ 
ter in our case does not have the terminating double- 
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quantum vortex observed in Ref. US]. Instead, there is 
one single-quantum vortex that reverses its circulation 
at the monopole, a scenario that has previously been 
shown to minimize the mean-field energy in the case of 
a Dirac monopole [37]. We confirm the presence of this 
single-quantum vortex in Fig. where it is visible as a 
line of suppressed spin density. We have verified that 
the phase winding along this vortex line reverses its sign 
near the origin where the magnetic field vanishes (data 
not shown). The orientation of the vortex depends on 
the particular realization of the applied noise. 

Figureshows the nematic vector and selected compo¬ 
nents of the spin vector during the decay of the isolated 
monopole. A ferromagnetic ring with a well-defined po¬ 
larization is clearly visible just after the creation ramp 
[see Fig.j^c) andj^d)] although it is so small that it was 
not distinguished within the finite experimental resolu¬ 
tion of Ref. [27] . This ring resides at the monopole core 
and retains its size during the temporal evolution. Since 
the nematic vector is observed to rotate by tt about the 
ferromagnetic core [see Fig. |^a)], the ring is identified 
as a half-quantum vortex ring, or Alice ring, discussed 
in Ref. [5]. We determine the radius of the Alice ring 
from the behavior of the nematic vector to be roughly 
0.2 jam which exceeds neither the spin healing length 
h/y/2m\c2\n{0) ^ 4 jam nor the density healing length 
h/y^2m\c2 + co|n(0) « 0.3 jam. Thus the structure es¬ 
sentially manifests itself a point defect. The subsequent 
decay of the polar phase destroys the Alice ring and even¬ 
tually the characterization of the condensate using the 
nematic vector becomes obscure. We therefore do not 
show the nematic vector but rather the local spin for 
long evolution times. Ultimately, the local spin aligns 
with the external magnetic field as shown in Fig. is) 
andj^h). 

Figure [^ shows the fraction of particles in the neu¬ 
tral state and the deviation of the order parameter from 
the initial isolated monopole state during the decay. 
The relative population of the neutral state is given by 
rinsit) = / dr|nns(r,t)'0(r,t)p, where A^(t) is the to¬ 

tal number of atoms and nns(i*,t) is a projection to the 
neutral state [13]. The deviation of the order parameter 
from the initial isolated monopole state at t = 0 is charac¬ 
terized by e{t) = 1 — 1/ dr'0*(r, t)'0(r, 0)\/ 

We observe that just after the creation ramp, roughly 
90% of the atoms reside in the neutral state, in agreement 
with the experimentally obtained value in Ref. [27] . The 
decay of the isolated monopole into the ferromagnetic 
phase is observed to change from an exponential-like be¬ 
havior into approximately linear decay with increasing 
magnetic field gradient. This observation suggests that 
a cascade of decay channels plays a significant role at 
strong field gradients. Due to the decreasing spatial over¬ 
lap between the resulting domains with increasing field 
gradient, the decay dynamics is slower the stronger the 


field gradient is. For a sufficiently strong gradient, the 
resulting domains are spatially well separated (see Sup¬ 
plemental Material). 



FIG. 4. Temporal evolution of the relative population of the 
neutral state, nns{t) (solid line), and deviation of the order 
parameter from the created isolated monopole at t = 0, e{t) 
(dashed line). In addition to the standard parameters given 
in the text, neutral state populations are also shown for (i) 
^ = 0 (O), (ii) r = 0 (-h) and (hi) C 2 = —A7Th^{a2 — ao)/(3m) 
(<C>) in Eq. M. Furthermore, we show nns(t) for cases, in 
which the quadrupole field gradient is suddenly changed at 
t = 0 from its standard value = 3.7 G/cm to 1.4 x (□), 
0.5 X 6q (x), and 0.3 x 6q (*). All curves are interpolated 
using cubic splines for enhanced visual appearance. 

Figure]^ also shows the results obtained for the three 
additional cases: (i) eliminating the quadratic Zee- 
man potential, (ii) eliminating three-body recombina¬ 
tion, and (iii) reversing the sign of the spin-spin interac¬ 
tion strength. None of these changes leads to a significant 
effect on the decay dynamics, indicating that the decay 
is not originating from these terms. We also studied the 
creation and decay of the isolated monopole with param¬ 
eters corresponding to ^^Na atoms and obtained qualita¬ 
tive agreement with the case of ^^Rb atoms (results not 
shown). 

Conclusions —Our numerical studies suggest that the 
isolated monopole structure observed in Ref. m con¬ 
tains a small Alice ring [5] . This vortex ring is destroyed 
by a subsequent dynamical phase transition into a ferro¬ 
magnetic order parameter supporting a Dirac monopole. 
Although the quadrupole field has been observed to stabi¬ 
lize the polar phase of a ^^Rb condensate if the field zero 
is well outside the condensate m, our simulations reveal 
that after the field zero is brought into the condensate, 
the polar phase decomposes on a time scale of 100 ms. 
We attribute this behavior to the spatially varying mag¬ 
netic field and the linear Zeeman interaction. Neither 
the spin-spin interactions, quadratic Zeeman effect, nor 
three-body recombination have a significant effect on the 
decay dynamics. However, the strength of the magnetic 
field gradient is shown to have a detrimental effect on 
the decay speed and characteristics. These studies set 
the stage for the detailed dynamics of topological point 
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defects in quantum fields. Finding ways to extend the 
lifetime of the defect further and thereafter to study the 
dynamics of multiple interacting point defects remain fu¬ 
ture challenges. 
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In Fig. SI, we show spin-contrast images of the condensate particle density for several waiting times after the 
creation ramp. The monopole conhguration in Fig. SI is created and evolved with a strong magnetic held gradient 
Rq = 11.1 G/cm which separates the SFSS and NS domains along z. 
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Figure SI: (a) Horizontally (b) and vertically integrated particle densities produced in the strong magnetic held 
gradient = 11.1 G/cm. All other parameters assume identical values to those specihed in the main text. The 
peak column density is hp = 2.7 x 10^^ cm“^ and the held of view is 22 x 22 (iro? in each panel. 




